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ABSTRACT 

The paper presents a method to recover exponential accuracy at all points (including at the 
discontinuities themselves), from the knowledge of an approximation to the interpolation polyno- 
mial (or trigonometrical polynomial). We show that if we are given the collocation point values (or 
an highly accurate approximation) at the Gauss or Gauss-Lobatto points, we can reconstruct an 
uniform exponentially convergent approximation to the function f{x) in any sub-interval of analyt- 
icity. The proof covers the cases of Fourier, Chebyshev, Legendre, and more general Gegenbauer 
collocation methods. 
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1 Introduction 


We continue our investigation of overcoming the Gibbs phenomenon, i.e., recovering pointwise 
exponential accuracy at all points including at the discontinuities themselves, from the knowledge 
of a spectral partial sum of a discontinuous but piecewise analytic function. In previous papers [3], 
[4], [5] and [6] we have considered the cases of Galerkin approximation, which means that, given 
an orthogonal basis over [—1, 1] with a weight function tu(x), we assume that the first N 

Galerkin coefficients 

1 

A = — / t v(x)f(x)<f>k(x)dx, 1 < k < N (1.1) 

Ck J - 1 

of a discontinuous but piecewise analytic function /(z), are known. Here c k are normalization 
constants. A procedure is given in our previous papers to recover uniform exponential accuracy, 
in any sub-interval of analyticity of f(x). The Fourier Galerkin case (<f>k(x) = e iknx , w(x) = 1) 
with a single discontinuity is analyzed in [3]; wave resolution issues are discussed in [4]; sub- 
interval reconstruction suitable for multiple discontinuities are analyzed in [5] for the Fourier case 
and the Legendre case (<f>k(x) = Lk(x), w(x) = 1), and in [6] for the general Gegenbauer case 
= C£(x), w(x) = (1 - x 2 y~ 2 ). 

The reconstruction procedure analyzed in these papers consists of two steps: 

1. Using the given spectral partial sum of the first N terms, to recover the first m ~ N 
Gegenbauer expansion coefficients, based on a sub-interval [a, b] C [—1,1] in which the 
function is assumed analytic, with exponential accuracy. This can be achieved for any 
L\ function, as long as we choose A in the weight function of Gegenbauer polynomials 
to be proportional to N . The error incurred at this stage is called the truncation error. 

2. For an analytic function in [a, 6], proving the exponential convergence of its Gegenbauer 
expansion, when the parameter A in the weight function is proportional to the number 
of terms retained in the expansion. The error at this stage is labeled the regularization 
error . 
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In this paper, we will consider the case of collocation. We assume that the point values f{xi), 
where Xi are the Gauss or Gauss- Lobatto points of the orthogonal basis {<j) k {x)}, of a discontinuous 
but piecewise analytic function f(x), are given. Alternatively, we assume that we have an accurate 
approximation to the interpolation polynomial based on these points. The objective is to recover 
exponentially accurate point values over any subinterval [a, 6] of analyticity of f{x). As before, 
we will separate the analysis of the error into two parts: truncation error and regularization error. 
Truncation error measures the difference between the exact Gegenbauer coefficients with A ~ A, 
and those obtained by using the collocation point values. This will be investigated in Section 3. 
The regularization error measures the difference between the Gegenbauer expansion using the first 
m ~ A Gegenbauer coefficients with A ~ A, and the function itself, in a sub-interval [a, 6], in which 
the function is assumed analytic. This error is estimated in [5] and we will simply quote the result 
in Section 4. The results are summarized in Theorem 4.2. Section 5 contains a numerical example 
to illustrate our results. In Section 2 we shall give some useful properties and estimates, including 
a useful approximation result. 

Throughout this paper, we will use A to denote a generic constant or at most a polynomial 
in the growing parameters, as will be indicated in the text. It may not be the same at different 
locations. 

2 Preliminaries 


We first collect some useful results about the Gegenbauer polynomials, to be used later. We rely 
heavily on the standardization in Bateman [1]. 

Definition 2.1 . The Gegenbauer polynomial C£(x), for A > 0, is defined by 

(1 - = G(A,n)^ [(1 - x a )" +A -i] (2.1) 


where G(A,n) is given by 


_ (~l) n r(A+ i)T(n + 2A) 
2 n n!r(2A)r(n + A + i) 


( 2 . 2 ) 
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for A > 0, by 


(7(0, n) 


. tirwi 

2 n-1 nr(ra + A) 


for A = 0 and n > 1, and by 

<3(0,0) = 1. 

for A = 0, n = 0. Notice that by this standardization, (7°(:r) is defined by (see [1]): 

C n( x ) = jC^(ar) = ^ T n (x ), n > 0; C’g(x) = 1, 

where T n (x ) are the Chebyshev polynomials. 


Under this definition we have, for A > 0, 


c£( i) = 


T(ra + 2A) 
n!T(2A) 


and 


|C»(*)I <C*(1), -1<x<1 

The Gegenbauer polynomials are orthogonal under their weight function (1 — x 2 ) A ~ 2 : 

j\ 1 - = 6 k , n h* 

where, for A > 0, 

We will need to use heavily the asymptotics of the Gegenbauer polynomials for large r 
For this we need the well-known Stirling’s formula: 

0 27r)h T+ 2e~ x < r(x + 1 ) < ( 2 tt )2 x x+ 2 e ~ x e^, x > 1 

Lemma 2.2 . There exists a constant A independent of A and n such that 

A2 


i-1 


■Cn(l) < h* < A- 


A2 


(n + A) '' nv ' / - "" - "(n + A) 
The proof follows from (2.9) and the Stirling’s formula (2.10) 


c n A (i) 


(2.3) 

(2.4) 

(2.5) 
□ 

( 2 . 6 ) 

(2.7) 

( 2 . 8 ) 

(2.9) 
and A. 

( 2 . 10 ) 

( 2 . 11 ) 
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□ 


Let us assume that the function f(x) is analytic in [— 1, 1] and satisfies the following 
Assumption 2.3. There exists a constant p > 1 and C(p) such that, for every k > 0, 

max 


*f (x) 

dx k 


k' 

< C(p W 
P * 


( 2 . 12 ) 

□ 


This is a standard assumption for analytic functions, p is the distance from [—1,1] to the nearest 
singularity of f(x) in the complex plane (see for example [7]). For such function f(x) we have the 
following 

Lemma 2.4 . If f(x) is an analytic function satisfying the Assumption 2.3, then for any integer m 
such that 0 < m < A — 5, we have 


max 

— 1 <X<1 


dx r 


0 -z 2 )M/(*)] 


<c W r(A + i)(2 + i)”. 


(2.13) 


Proof: 


dx T 


(.-* 2 > A -i/(*>]| = (&*>) 

ml (m - k)\ 2 fc fc!r(2(A-fc))r(A+l) 


£ k '( m - *) ! P m ~ k r (a - k + i) r( 2 A - k) 

< r(n^ m! 

- p* 


] )\ 

(l-x^- k -i\c x k - k (x)\ 


2 k k\ T ( 2 (A — k)) r | 

( A+ 2) 

1 r(2A- k) 

r | 

[a - k + \ j 

| r(2A - k) 

Ar! T (2(A - k)) 


< C(p) 2 m r(A + -)x: 


m! 




< 


2/ *=o fc! r (a - * + i) V 2 a> 
cW2"r(A + i)£^^(±) 
c W r(A + I )( 2 + J 


m—k 


where we have used the analyticity assumption (2.12), and the formulas (2.1)-(2.2), for the second 
inequality. We have also used (2. 6)-(2. 7), and the fact that (1 — x 2 )^~ k ~ 2 < 1 because A — k— | > 0, 
for the third inequality. In the fourth inequality we have used the fact that Y ^A — k + 5) £ (m-k)! 
because A — k + £ > m — k + 1. 
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Finally, we would need to quote some approximation results for the interpolation polynomials. 

If the point values /(x,) of the function /(x) are known at the N Gauss or Gauss-Lobatto points 

{x t } of the basis functions {<^(z)}, we will denote the unique interpolation polynomial by 

N 

InHz) = 53 Mk(x) (2.14) 

k - 1 

which satisfies 

iNf(xi) = /(*.')» i=l,2,...,N (2.15) 

We then have the following error estimates for the interpolation polynomial Ijyf(x): 

Lemma 2.5 . If /(x) has m continuous derivatives in [-1,1], then the interpolation polynomial 
In f(x) defined by (2.14), (2.15), where the basis functions {<^(x)} are either trigonometric poly- 
nomials e lk ' KX or the Gegenbauer polynomials with ft > —\ (which includes, for example, the 
Chebyshev case fi = 0 and the Legendre case fi = ^), satisfies the following estimate: 

ll/-/*/ll«<}^ll/ ( " l lk-. ( 2 - 16 ) 

where the weighted L 2 norm is defined by 

\\f\\h = [ w(x)\f(x)\ 2 dx (2.17) 

and A is a constant independent of N and m. 

Proof: The proof of this Lemma can be found in [8] for the trigonometric polynomial case, and 
in [2] (Theorem 6.6.1, first line of formula (6.6.3), for which the only restriction about the weight 
corresponds to our fi > — ^) for the Gegenbauer polynomial case. Notice that in both situations 
the constant A is actually independent of m, although the authors did not explicitly point this out. 

□ 


3 Truncation Error 

Consider an function /(x) defined in [—1, 1]. We assume that the point values f{x%) at the N 
Gauss or Gauss-Lobatto points (i.e., the zeros of <f>w+ i(x), or the zeros of j^4>n{x) plus the two 
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boundary points ±1), are given. Here, the orthogonal basis {</>*($)} is either the trigonometric 
polynomials e lk * x , or the Gegenbauer polynomials C^{x) with p > — 

We are interested in approximating the Gegenbauer expansion of f(x ), with A ~ A r , based on 
a sub-interval [a, b] C [—1, 1] of analyticity. We start by introducing the local variable 

Definition 3.1 . The local variable £ is defined by 

£ = f(z) = > X = x(f) = ff + tf (3.1) 


where 


e 


b — a 

2 ’ 


Thus when a < x < b y — 1 < £ < 1. 



(3.2) 


□ 


We assume that the function 


/(x) is analytic in this sub-interval [a, 6] satisfying the Assumption 


2.3: 

max 

a < j ;<6 

for any k > 0, with some p > 1. 

The function /(a;) has a Gegenbauer expansion in the sub-interval [a, 6], with A ~ TV. With 
6 and £ defined in (3.1)-(3.2), we have 

oo 

f(tt + 6) = Y / f<(l)CKO, -1<£<1 (3.4) 

1=0 

where the Gegenbauer coefficients /*(/) are defined by 

/?(<) = A/j 1 - « 2 ) M C, A (f)/M + S)<K (3.5) 

Our goal in this section is to find a good approximation to the first m ~ N Gegenbauer 
coefficients f^(l) in (3.5). 

Based on the known point values /(#;), we define not the usual interpolation polynomial Iwf(x) 
in (2.14)-(2.15), but 

g N (x) = I N (a */)(x), (3.6) 


d^£ 

dx k 


(x) 


- C ^~pk 


(3.3) 
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where the function a(x) is defined by 


a(x) 


= f (i-t(z) 2 ) 

\ 


2'iA-i 


a < x < b, 
otherwise. 


(3.7) 


with £(x) defined by (3.1). This works for the case dx < oc, where w(a:) is the weight 

function of the basis Most commonly used bases belong to this class, for example the Fourier 

case where w(x) = 1 and the Gegenbauer case where ie(x) = (1 — x 2 )^~^ , when fi < Otherwise, 
if the basis {<t>k} consists of the Gegenbauer polynomials C^(x) with /x > we should choose 


a(x) 


(1 — £(x) 2 ) a 2/(1 — x 2 )% a < x < 6, 

0, otherwise . 


(3.8) 


Intuitively, the function gi\r(x) has about A ~ N continuous derivatives, hence would produce 
nice error estimates. 

Our candidate for approximating the Gegenbauer coefficients / A (/) in (3.5) is: 

i r x 


9<(l)= ^ J^gN(*t + 6)CH Orff 

if (3.7) is used to define a(x), and is 

gHn = A- f{ i - 

if (3.8) is used to define a(x). 

How well do # A (/) approximate / A (/)? To answer this question we define 
Definition 3.2 . The truncation error is defined by 

m 

TE{ A, m, N, c) = max | £(/ A (/) - $ € A (/))C?(OI 
S * S /=0 

where / A (/) are defined by (3.5) and g^(l) are defined by (3.9). 


(3.9) 


(3.10) 


(3.11) 


The truncation error is the measure of the distance between the true Gegenbauer expansion in 
the interval [a, 6] and its approximation based on the collocation point values of f(x) in [—1, 1]. 
Using the result of the previous section, we have the following estimate for the truncation error: 
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Lemma 3.3 . The truncation error can be estimated by 


TE(\,m, N,() < A 


m(m + A)f ^ T(m + 2A) p ^ 


v / Am!r(2A)^ A -t-f 

Proof: We shall just prove the case when (3.7) and (3.9) are used. For this case we have: 

m 

TE(\,m,N,<) < El/"(0-ff f A (0IQ A (l) 

/=0 

< p \j\ [(1 - £ 2 ) A “> /K + 6) - 9n(<S + 0] CHOdi 

m fc*/ A (l)) 2 i r * 

< 7T\ ~ 7 / Hx)f{x)-g N {x)\dx 

7^0 h i ( J - 1 

^ E fl/ (y_ i W(a:)|Q( a :)/( a: )-g j v(x)| 2 dx) Q ^ ^dx 

. f>( C "(l)) 2 A 


(3.12) 


< 


< A 


Ml 

br £ 


i r H)K)‘l 


Nl x ~tt t'o 

m(m + A)T (a + ^ T(m + 2A) j-j ^ 


y/X ml r(2A)7V A -2 

where we have used (2.7) for the first inequality; the definitions of f*(l) and <7^(/), in (3.5) and (3.9), 
for the second inequality; the definition of a(x) in (3.7) for the third inequality; Cauchy-Schwartz 
inequality for the fourth inequality; the approximation result (2.16) for the fifth inequality, with 
[A - |] denoting the greatest integer < A - (2.13) for the sixth inequality; and (2.6) and (2.11) 

for the last inequality. 

The case when (3.8) and (3.10) are used instead of (3.7) and (3.9) can be proven in a similar 
fashion. We only need to notice that the function 2 ? > s always analytic in a < x < 6, 

_ £i 

hence can be combined with f(x) in the analysis. This accounts for the extra factor N 2 in the 
denominator. 


□ 
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We are now ready for the main theorem of this section: 

Theorem 3.4. Let the truncation error be defined in (3.11). Let A = acN and m = j3(N with 
0 < a,0 < 1, then 


TE(a(N,/3(N,N,e) < A 


'(/3 + 2af +2a (\ 1 \ a V 


e a a a (3P 


(K)') 


tN 


where A grows at most as JV 2 . In particular, if a = 3 < —r ?e , v , then 

27 ( 1 +^) 

T E(aeN, atN, N, e) < Aq< N 


where 



Proof: The theorem follows from Lemma 3.3 and the Stirling’s formula (2.10). 


(3.13) 


(3.14) 


(3.15) 


□ 

We now make two remarks: 

Remark 3.5. Unlike in the Galerkin case [3], [5] and [6], where for the purpose of truncation 
error being exponentially small we only need to assume that the function f(x) is in L\ here for 
the collocation case we need to assume that the function f(x ) is actually analytic in the relevant 
sub-interval [a, b]. 


□ 

Remark 3.6 . We would like to point out that the collocation procedure (3.6), where we collocate 
not the function f(x) but its product with a weight (l-x 2 ) A ~?, is somehow nonstandard. Another, 
probably more natural, procedure of collocation, is to interpolate f(x) with Ijvf(x) and then define 
the approximate Gegenbauer coefficients <? A (/) by 

9 ? (/) = -nfj 1 - (3.16) 

instead of (3.9). Unfortunately, we have been unable to prove that the truncation error is expo- 
nentially small for (3.16). In fact it can be shown that for some cases the aliasing error reduces 
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the accuracy. However numerical experiments in Section 5 seem to indicate that both procedures 
produce similarly good results. 


□ 


4 Regularization Error and the Main Theorem 


The second part of the error, which is called the regularization error and is caused by using a 
finite Gegenbauer expansion based on a sub-interval [a, 6 ] C [—1, 1], to approximate a function f(x) 
which is assumed analytic in this sub-interval, has been studied in [5]. We will thus just quote the 
result. 

We again assume that f(x) is an analytic function on [a, 6 ] satisfying the Assumption 2.3, see 
(3.3). 

Let us consider the Gegenbauer partial sum of the first m terms for the function /(c£ + S): 


rn 

/£■({) = E/*(')ef(0 


/= 0 


(4.1) 


with ( and 6 defined by (3.1) and (3.2), and the Gegenbauer coefficients /*(/), based on [a,b], 
defined by (3.5). 

The regularization error in the maximum norm is defined by: 


RE(\, m, e) 


max 

— 1<£<1 


J=0 


(4.2) 


We have the following result for the estimation of the regularization error, when A ~ m [5]: 

Theorem 4.1 . Assume A = 7 m where 7 is a positive constant. If f(x ) is analytic in [a, b] C [-1,1] 
satisfying the Assumption 2.3 in (3.3), then the regularization error defined in (4.2) can be bounded 

by 


RE^m, m,e) < Aq m 


(4.3) 


where q is given by 


_ c(l + 2 7 ) 1 + 2 t 

q ~ / 92 1+2 ^(1 + 7 ) 1+ ' 1 ' 


(4.4) 
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which is always less than 1. In particular, if 7 = 1 and m = /3eN where /? is a positive constant, 
then 


with 


RE(fieN, (3tN,c) < Aq 


.(**)' 

\32pJ 


:N 


(4.5) 


(4.6) 


□ 


We can now combine the estimates for truncation errors and regularization errors to obtain the 
following main theorem of this paper: 

Theorem 4«2 . (Removal of the Gibbs Phenomenon for the collocation case). 

Consider a L\ function f(x) on [—1, 1], which is analytic in a sub-interval [a, b] C [—1,1] and 
satisfies the Assumption 2.3 in (3.3). Assume that the point values f(xj) at the N Gauss or Gauss- 
Lobatto points of {<fik( x )}, are given. Here {<^(2;)} is either the trigonometric polynomials e xkirx or 
the Gegenbauer polynomials C£(x) for fi > — Let #*(/), 0 < / < m be the Gegenbauer expansion 
coefficients, defined in (3.9) or (3.10), based on the sub-interval [a, 6], of the collocation polynomial 
9n{x) defined in (3.6). Then for A = m - /3eN with / 3 < 2e , T , we have 


27 0+£)’ 


max 
-1 <£<1 


171 

+ *)- £$*(/)<?(£) < A (qf + q' R ") 
1=0 


(4.7) 


where 


<IT = 


m f 1 + i 

2e \ + 2pJ\ 


10 


< 1, 


QR 


m* < . 


V32 p) 


r 5±i £ 


and A grows at most as A 2 
Proof: Just combine the results of Theorems 3.4 and 4.1. 


We remark tnat no attempt has been made to optimize the parameters. 


□ 
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5 Numerical Results 


In this section we give a numerical example to illustrate our result. We will test both the Fourier 

and the Chebyshev cases. 

Example 5.1 . We take the following function 

f(x) = / sin(cos(x)), if a < x < b (5.1) 

J ' ' I 0, otherwise 

and try to recover pointwise values of f(x) for a < x < 6. This function is also used in [6]. In order 
to see the effect of relative locations of the sub-interval, for the Chebyshev case we will consider 
two situations [a, 6] = [-0.5, 0.5] and [a, b] = [0, 1]. For the Fourier case we will only consider the 
case [a, 6] = [-0.5, 0.5] because of periodicity. 

We first consider the Fourier case. We assume that the point values /(x;), for the 2 N + 1 equally 


spaced points: 


2i 


x, = 


-N <i< N, 


(5.2) 


2N + 1’ 

are given. Two collocation procedures are tested. Procedure A, which is the procedure proven in 
Section 3 to be exponentially accurate, is to interpolate a(x)/(x) as in (3.6) and (3.7), and then 
build the approximating Gegenbauer coefficients g^(l) by (3.9), finally use 


m 

Sm(l) = 5Z 9t( l ) C HO ( 5 ‘ 3 ) 

/=0 

as our approximation to /(x) in [a, 6]. We also test the more natural collocation Procedure B , 
for which we interpolate /(x) as in (2.14)-(2.15), and then build the approximating Gegenbauer 
coefficients g*(l) by using (3.16), finally using (5.3) to approximate /(x) in [a, b]. 

For the Fourier method, since the collocation points are equally spaced, we also have the choice 


to use all the collocation points in [-1,1] to build the approximation in [a, 6] (referred to as global- 
local)] or use just the collocation points inside [a, 6], with suitable scaling to expand [a,h] to [—1, 1] 
(referred to as local- local). 

We shall use the parameters 


m = 0.2(N, A = 0.4efV 


(5.4) 
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where 2N + 1 is the total number of collocation points in the full interval [-1, 1]. This choice is also 
used in [6] (with the right scaling). We will use (5.4) for all the test cases of Fourier collocation. 
No attempt has been made to optimize these parameters for each individual case. 

In Fig. 1 we show the errors, in a logarithm scale, for N = 10, 20, 40 and 80, for the Procedure 
A ( global-local at left, local-local at right). In Fig. 2 we show similar errors for the Procedure 
5. We can clearly see that both Procedure A and Procedure B produce similarly good results. 
Also, it seems that global-local and local-local are producing similar errors for the same choice of 
A and m. Considering that the effective N (the number of collocation points) is smaller for the 
local-local case, hence round off errors should also be smaller, it seems that local-local should be 
used if possible (see also the discussions in [9]). However, there are situations (for example when 
we are solving a PDE and 7yv/ is obtained from the time evolution) where the point values f(x t ) 
may not be accurate, however the interpolation polynomial I N f (which may interpolates f(x) at 
some other points than is accurate. In such situations global-local is more appropriate. 


Fourier, Global- Local, Collocate Pweight 



Fourier, local-local, collocate f*weight 

to 1 r 
'o’ r 

I0 4 r 

10 4 - 

10' r 
10 4 r 
10* r 
io- r 

to" r 
10 " r 
io“ r 

ur M »■■ ■■! .... i .... l .... t .... i .... i 

-0.4 -0.2 0.0 0.2 0.4 00 



Fig. 1: Errors in log scale, f(x) defined by (5.1). [a, 6] = [—0.5, 0.5]. Fourier collocation with 
Procedure A, namely interpolating f(x) * a(x) with a(x) defined by (3.7). A = 0.4eiV and r n = 
0.2(N, N = 10,20,40,80. Left: global-local ; Right: local-local. 
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Fig. 2: Errors in log scale, /(x) defined by (5.1). — [—0.5, 0.5]. Fourier collocation with 

Procedure J 9, namely interpolating f(x) only. A = QAeN and m = 0.2 eN, N = 10,20,40,80. Left. 
global-local ; Right: local-local. 


Next we consider the Chebyshev case. We assume that the point values /(#t)> the N + 1 
Gauss- Lobatto points 

= cos ^ ? 0 < i < N, (5.5) 

are given. We again test both Procedure A and Procedure 5, now for the middle interval [a, 6] = 
[-0.5, 0.5] as well as for the boundary interval [a, 6] = [0, 1]. The choice of A and m are the same 

as in [6]: 

m = 0.UN, A = 0.2 eN (5.6) 

This is also the same as (5.4) for the Fourier case, considering that now the number of collocation 
points is N + 1 instead of 2 N + 1. 

In Fig. 3 we show the errors, in a logarithm scale, for N - 20, 40, 80 and 160, for the 
Procedure A ([a, 6] = [-0.5, 0.5] at left, [a, b] = [0,1] at right). In Fig. 4 we show similar errors 
for the Procedure B. We can again see that both procedure A and Procedure B produce similarly 
good results. Also, these errors are comparable with those obtained in [6] where Galerkin instead 
of collocation is used. 
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Fig. 3: Errors in log scale, f(x) defined by (5.1). Chebyshev collocation with Procedure A, namely 
interpolating /(*)*«(*) with a(x) defined by (3.7). A = 0.2 eN and m = O.lcTV, N = 20,40,80,160. 
Left: [a, 6] = [-0.5, 0.5]; Right: [cf,fc] = [0, 1], 




Fig. 4: Errors in log scale, f(x) defined by (5.1). Chebyshev collocation with Procedure B, namely 
interpolating /( i) only. A = 0.2(N and m = O.lcJV, N = 20,40,80,160. Left: [a, 6] = [-0.5, 0.5]; 
Right: [a, 6] = [0, 1]. 


6 Concluding Remarks 

We have proven the exponential convergence in the maximum norm, of a reconstruction procedure 
using Gegenbauer series based on Cf( x) with large A, for any Li function in any sub-interval [a,h] 
in which the function is analytic, if we are given the collocation point values of the function at the 
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Gauss or Gauss-Lobatto points. The proof covers Fourier, Chebyshev, Legendre and more general 
Gegenbauer collocation methods. Numerical examples are also given. 
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